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Abstract 

We study an integro-differential equation which generalizes the periodic interme- 
£NJ ■ diate long wave (ILW) equation. The kernel of the singular integral involved is an 

elliptic function written as a second order difference of the Weierstrass ^-function. 
Using Sato's formulation, we show the integrability and construct some special so- 
lutions. An elliptic solution is also obtained. We present a conjecture based on a 
Poisson structure that it gives an alternative description of this integrable hierarchy. 
We note that this Poisson algebra in turn is related to a quantum algebra related 
with the family of Macdonald difference operators. 
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1 Introduction 



In this paper, we consider an integrable differential equation with a singular integral term 
associated with a doubly periodic function. We often classify the known integrable equa- 
^-j- ■ tions with singular integrals according to the periodicity of the kernel functions. Namely, 
the Benjamin-Ono equation [TJ [2] corresponds to the case with no period (rational func- 
tion) because it has the Hilbert transformation, and the intermediate long wave (ILW) 
equation 0, H] corresponds to the singly periodic case (trigonometric function) as hyper- 



bolic cotangent is involved as the kernel. The doubly periodic case was introduced in [5] 
as a periodic version of ILW equation. 

We aim at constructing an integrable equation which recovers all these equations as 
special limits, and which also relates to the theory of the Macdonald polynomials [6] 
(see [3 [8] also). We propose that the kernel of our singular integral is an elliptic function 
having simple poles at three points 7,0, —7 in the fundamental parallelogram with residues 
1, —2, 1 respectively and holomorphic elsewhere. 

Let ui, u 2 and 7 be complex numbers such that the ratio 5 = U2/0J1 satisfies Im(<5) > 0, 
and < Im(7) < Im(<5). Let x and t be real independent variables, and let rj(x,t) be an 
analytic function satisfying the periodicity condition rj(x + l,t) = r](x,t). We consider 
the integro-differential equation 

d iui f 1 ^ 2 

- v (x,t)= V (x,t)-^ { (A 1 ()(2uj 1 (y-x))- V (y,t)dy, (1) 

Ot 7T J_ 1/2 
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where ((x) = 2ui, 2o; 2 ) denotes the Weierstrass C-function [9] (see f|96|) in Appendix 
A), the discrete Laplacian A 7 is defined by (A ry f)(x) = f(x — 7) — 2f(x) + f(x + 7), and 
the integral J means the Cauchy principal value. 

Our purpose in this paper is to study (CD) from several viewpoints. First we use the 
standard method which transforms ([1]) into a difference-differential form [5j [TUJ, [TTj, [12]. 
Then we follow the method developed in [131 EE] to utilize the Sato theory to construct an 
integrable hierarchy which includes ([!]). Then we study the system of integrals of motion 
associated with (jTJ) in terms of the Sato theory. 

We remark the following. In [TQj, [T2], the conserved densities for the (periodic) ILW 
was studied by using the Backlund transformation. However, we do not know a Backlund 
transformation for ([Tj) at present. It is an open question to find it and compare the 
approach given in this paper with those classical analysis of conserved quantities. 

In the papers [T5J rTB] , a Hamiltonian approach to the ILW equation was pushed forward 
by using the Gel'fand-Dikij brackets and the bi-Hamiltonian structure. We will develop 
a Hamiltonian description of ([!]) in the same spirit as theirs. Our situation, however, 
might be a little tangled in the following sense. In one hand, it has a direct connection 
with the Macdonald difference operators [6], or to be more precise, to its elliptic analogue 
defined through the algebra of Feigin and Odesskii [T7| (see [7] and [8]). On the other 
hand, we also attempt to connect the Hamiltonian approach with the Lax formulation of 
Sato. It is a future problem to understand analogues of the Gel'fand-Dikij brackets and 
the bi-Hamiltonian structure for (pQ). 

Finally, we make an important comment that the Poisson algebra in this paper has a 
deep connection with the one found in [18J. The difference analogue of iV-th KdV studied 
by Frenkel has two parameters q and N. It can be found that if we set q = e 2mj and 
impose the condition 5 = N^y, we almost recovers Frenkel's Poisson algebra, but missing 
the delta function terms which typically appears in the deformed W-algebras (see [19] ). 

This article is organized as follows. In Section 2, we rewrite ([1]) in a form of differential- 
difference equation. Then we show the ordinary ILW equation with periodicity can be 
obtained from ([1]) in the limit 7^0. In Section 3, by using the standard Sato theory 
we present an integrable hierarchy which contains ([T]) in the lowest order. We study 
the structure of the integrals of motion in some detail in this setting. Finally, Section 
4 is devoted to an alternative description based on a Poisson structure derived from a 
quantum mechanical integrable model associated with the Macdonald theory, from which 
we recover the same equation ([I]), and presumably all the equations given in the hierarchy. 



2 



2 Differential-difference form 



2.1 integral operator T 

Let T be the integral transformation denned by 

(Tf)(x) = ^ / 1/2 (A 7 C)(2^i(y - x)) ■ f{y)dy, (2) 

Tf 7-1/2 

then, (TjQ) can be written as fj = rj (Trj). 

Decomposing r](x, t) by the Plemelj formalism [20], one obtains a differential-difference 
equation and a bilinear equation from ([1]). Here and hereafter, we set p = e 2jTlS , q = e 2nn 
for simplicity, and let D be a domain in the complex z-plane containing the infinite strip 
< Im(z) < lm(5). 



Lemma 2.1 For any nonzero integer m we have 



f ] — a m \('\ — rt m n~ m \ 



V 



and T(l) = 0. 



Proof. Use ( )97|) in Appendix A. □ 

Corollary 2.2 Suppose that an analytic function f(z) on D satisfies the periodicity f(x+ 
1) = f(x), then we have 

T(f(x) -f( x + 8)) = f(x) - f{x + 7 ) - fix + 5 - 7 ) + fix + 5). (4) 

Define difference operators T and S acting on the variable x by Tfix) = fix + 5), 
Sfix) = fix + 7). For simplicity of display, we also write / = Tfix) = f(x + S),f = 
Sfix) = fix + 7), / = S' _1 /(x) = fix — 7) and so on. Then we can write (j3J) as 

T(/-/) = (l-5)(l-T5- 1 )/, (5) 
for example. Setting g = (1 — T)/, this 'formally' can be expressed as 

Proposition 2.3 Suppose that w{z) is holomorphic on D and satisfies the periodicity 

1 /2 

wix + 1) = wix). Set rjix) = — w(x + 5) + r] , where r] = f_y 2 r]ix)dx denotes the 
zero Fourier component. Then we can recast (T7]j into the difference equation 

d _ 

— (w — w + 770) = (w — w + 770) (w — w — w + w) . (7) 
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Remark 2.4 Note that we have dr]o/dt = from the assumption on w(x). 
Proposition 2.5 Let e and r/ be constants. Assume that r satisfies the bilinear equation 

D t r ■ r = eff - t]qTT, (8) 

Q T 

where D t denotes the Hirota derivative defined by D t f-g = (dtf)g-f(dtg). Setw = . 

r 

Then w satisfies the difference equation 
Proof. From ([8]) we have 

W - W + 7/o = h + TjQ = £ — , 

T T TT 



hence 



on/ - \ <9t r <^ r ®tZ d t r . _ ^ 

Ot logiw — w + 1]o) = h + — — = (w — w — w + w). 

T T T T 



□ 



2.2 some special solutions 

We give some examples of special solutions to ([I]), ([7j), or (jSJ). 

Proposition 2.6 Let n be a positive integer. Let Ci, • • • ,c n be n complex parameters, 
ki, ■ ■ • k n be n integers which are all distinct and nonzero. Set 

t = det(fi >m )i<i >m < n , (9) 
fi,m = + cm™- 1 exp(2mkix + (m - A,)t), 

I _ e 27ri(<5-7)fc ; J _ e -27ri(5-7)fc; 

^ ~~ ~ £ I Z e 27ri<Sfc ; ' ^ ~~ — £ x — e ~ 27TiSkl ' 

T/ien t/iis tau satisfies the bilinear equation written for rjo — e. 

The proof of Proposition 12.61 will be given in Section 13.21 
Remark 2.7 We have 



Now we turn to the case of an elliptic solution. Let A be a complex number satisfying 
Im(A) > and define 

A) = _ i ' e 7riA / 4 + 7riu V^^_2^m e 7riAm(m+l)+27rim« ^-q 
■mGZ 



Proposition 2.8 Let e and A be a complex parameters satisfying Im(A) > lm(S), and 
let k be a nonzero integer. Set 

t = -di(kx + u(k)t, A), (12) 
^(fc 7 ,A)^(^-7),A) 

u{k) = ~ £ ' (13) 

and 

= »(" U«7-«).A) " ^A) J ■ (14) 

Then r satisfies (Ej) and 

V = e= 1 , (15) 

T T 

is a special solution to the integro- differential equation (QJ). 

It is straightforward to check this by using addition formulas, so we omit the proof. 

Remark 2.9 Note that, the A (not the 5 of the period of the equation (QJjj gives the 
period of the solution. 

2.3 limit to periodic ILW 

Now we consider the limit 7 — > 0, and derive the periodic ILW equatin studied in [5]. Let 
T be the integral operator 



(Tf)(x)= f / \((2u Jl (y-x))-2((u; 1 )(y-x)}-f(y)dy, (16) 

TT 7-1/2 



then we have 



1 1 +p m 



r(e 2„ } = (m ^ ), (17) 

and T(l) = 0. Hence we have the formal expression 

r s = -\\^ 9 . (18) 

Lemma 2.10 We have the expansion ofT in 7 as 

T = -7^ + 7 2 T9 2 + 0( 7 3 ). (19) 

To have the periodic ILW equation, we need to cancel the leading term — ^fd x . To this 
end we first perform a Galilean transformation on f](x,t) and assume the expansion as 
e ctdx r] = e + ju(x, t) + 0(7 2 ), where c and e are constants. Rescaling t as t — > 7~ 2 t and 
setting c = £7 + cry 2 + • • • , we have the integro-differential equation of ILW-type for u 

u t = au x -uu x + eT{u xx ). (20) 
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3 Lax formalism 



3.1 Sato theory 

We introduce an infinite set of independent variables x, r and t = (t\ = t, t 2 , t 3 , ■ ■ ■ ). 
Define difference operators S, R and T by Sf(x, r, t) = f(x + j, r, t), Rf(x, r, t) = f(x, r + 
6, t), Tf(x, r, t) = f(x + 5, r, t). For simplicity, we write Sf = f, Rf = f, Tf = f and so 
on. 

We formulate a version of Sato theory based on the papers [13] and [H]. In what 
follows, we work with a space of operators expressed as formal series in S~ x . Let W = 
1 + WiS^ 1 + w 2 S~ 2 + ■ ■ ■ be the Sato- Wilson operator. Let 0(A) and k(X) be Laurent 
series in A as 

<KA) = A + ^0 m A-\ /<A) = A + 5^ +1 A-*, 0i,/%eC, (21) 

j>0 i>0 

and define p(A) = </> -1 (A) by the condition 

p(A) = A + ft+i A"*, 0(p(A)) = A. (22) 

We impose the following set of evolution equations on W: 
dW 

+ jy </>(£)■> = = (W / 0(5)W" 1 )+, (23) 



at, 1 "> 

Wk{3) = CW, C = (Wk{S)W~ 1 ) + , (24) 
where we have used the standard notation that we write (A)+ = Yl°j=k a j^~'' f° r an y 

Remark 3.1 The equality p^] ) means the truncation of the operator W k{S)W~ 1 = 
(Wk(S)W~ 1 ) + , which involves an infinite system of equations about integrals of motion, 
and is very much important for our task. Namely, the K m 's plays the role of integrals. On 
the other hand, the parameters <p m 's are introduced just to have some linear combination 
of the time variables ti ? s. See Section EOl and Section\3J h . 



Example of the time evolution (1231) : We have 

fc-i 

d tl w k = w k (wi - Wi) - (w k+1 - w k+1 ) - <p k+ i - 22 4>j+i w k-j, (25) 

3=1 

etc. 
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Proposition 3.2 The compatibility among the evolutions in the t directions and the r 
directions can be described by the Zakharov-Shabat equations 

§|-|H B ^l = > < 26 > 

BC ~ 

— — h CBi - B % C = 0. (27) 
dti y ' 

Introduce the wave function \l/ defined by 

r, t) = Wp(\) xh K(X) r l s e^ x \ (28) 

oo 

e(t,A) = 5^t i A*. (29) 

i=l 

Proposition 3.3 T/ie u>ai>e function satisfies 

(30) 

* = Ctf . (31) 
The operators B\ and C can be written explicitly as 

Bx — S + w\ — w\ + 0i, C = 5" + W\ — W\ + K\. (32) 
Write w = W\ for simplicity, then from Zakharov-Shabat equation (I27P we have 



— (w — tu + «i) + (iu — w — w + w)(w — w + k{) = 0. (33) 

Definition 3.4 We call the condition on W 

W(x,r + 5,t) = W(x + j-5,r,t), (34) 
(W = T~ l SWTS~ x or w k = Wk) the reduction condition. 

Proposition 3.5 Assume that W satisfies the reduction condition (34\ ), and define the 
Lax operator by L = TS^C = WT S^ 1 k{S)W^ 1 . Setting r]o — K\, r\ = w — w + r] , we 
have L = T + rjTS~ l , and w satisfies the difference equation (7|). Hence r\ satisfies (CP 
under the condition that w is holomorphic on D and periodic w(x + 1) + w(x). 

Proposition 3.6 Assuming the reduction condition (34\ ), we have 

Lm = p(A) (5 - 7)/7 /t(A)*. (36) 



7 



3.2 Casorati determinant 

As an application of the Sato theory we have studied in the previous subsection, we 
construct a special solution of (JT]) in terms of a Casorati determinant. In this subsection 
we restrict ourself to the simplset possible situation 0(A) = p(A) = A and k(A) = A + e 
and assume the truncation of the Sato- Wilson operator as 

W = 1 + w^- 1 + w 2 S~ 2 + ■■■ + w n S~ n . (37) 

We note that the case e = necessarily gives us the trivial situation rj = 0. For simplicity 
we sometimes denote S k f(x) = p k \x) etc. 

Consider the linear system WS n f = 0. Let fi, f 2 , • • ■ , f n be the basis of the linear 
system satisfying the dispersion relation 

wr s ' f " (38) 

Rfi = (S + £)/,. (39) 
As usual, we explicitly give the basis as 




by introducing the set of parameters {Aj,//j, Cj|l < i < n}. 

Proposition 3.7 The operator W uniquely characterized by the linear system WS n fi 
(i = 1, • • • ,n) satisfies the evolution equations l[2!^) ,(24\) written for 0(A) = A,p(A) = 
A, k(X) = A + s. 

For any sequence of integers ai,a 2 , ■•■ ,a n , we denote the correspondong Casorati 
determinant by the symbol \ai,a 2 , ■ ■ ■ ,a n \ = det(// aj ' ) )i<i i j< n . Setting t = \n — l,n — 
2, • • • , 0|, the Wj's are written as 

w k = {-l) k \n,n- I,--- ,fc,-.. ,0|/r, (41) 

where the symbol k means that the letter k is eliminated. Introduce differential operators 
Pk by 

OO / OO \ 

J2 Pk y k = expl-J2-dt m y m )- (42) 

fc=0 V m=l / 
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Proposition 3.8 We have Wk = {Pk T )/ T , namely 

f;«>(»,r,t^= T(3! - r ; t -[f 1 > , (43) 

where we have used the standard notation [A] = (A/1, A 2 /2, A 3 /3, ■ ■ ■ ). 

Note that we especially have Wi = —(d tl r)/T, w n = (— l) n (SY)/r and w n+ i = w n+2 = 
• • • = 0. Hence we have 

B l = S + d -^, C = S + e^. (44) 



Proposition 3.9 If the condition 

Xi + e -(\ j M^\ (Jm 



fij + e 

is satisfied, the reduction condition (34\ ) holds. Such a pair is parametrized by via 
kj as 

I _ g27ri(5-7)fcj -y _ g — 2m(S—y)kj 

\ = ~ e i _ e 2niSkj ' ^3 = ~ e l _ e -2-KiSkj " (46) 

Assuming the reduction condition (13"4"|) and setting L = W(T + eTS'^W" 1 = T + 
rjTS -1 , we have the representation of the rj which satisfies ([T|) in terms of the r in two 

ways as 

V = \~ — + e = e^—. (47) 

T T T T 

This means that the r satisfies the bilinear equation (18"1). Thus we have the special solution 
presented in Proposition 12. 61 

The wave function \I/ associated with this solution reads 

* = r(x ' r ; t-[ f 1]) A^(A + ey'SW, (48) 
r(x,r,t) 



which satisfies (under the condition ([3] 

= A (5 ~ 7)/7 (A + e)*. (49) 
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3.3 conservation laws 

Throughout this subsection, we assume that the reduction condition (I3"lj) is satisfied. 
Introduce the Fourier expansions 

rj(x) = V- n e 2mnx , w k (x) = £ w k ,. n e 2 ™ nx , (50) 

nGZ ngZ 
rl/2 f l/2 

Vn = / dxr](x)e 2mnx } w k , n = / dxw k (x)e 2 ™ nx . (51) 

7-1/2 J-l/2 

Conjecture 3.10 Assume the evolution equation in the r direction [Efy and the reduction 
condition (34\ ). The K k 's are given as certain degree k expressions in {i] n }. In other words, 
r](x) satisfies infinitely many constraint conditions for given {n k }. 

Example: 

K t = T] , (52) 



'a _— n 



Ep q 
Y—^VnV-n, (53) 



Ep y p q sr^ p q ir .., 

T~n^ y—;VmV-m + nV-n ~ ^ Q _ ^ OVnV-n , (54) 

etc. 

We present some explicit calculations to explain what is meant by Conjecture 13. 10L 
Lemma 3.11 It follows from L = T + r]TS~ l = WTS^^W' 1 that 

fc-i 

nm.k-i = w k - w k + ^2 Kj+iWk-j-i- (55) 

j=0 

First we consider the case k — 1 in (151)1) . namely r] = Wi — w\ + K\. From this we have 
the relations among the Fourier modes of W\,T] and K\ as 

wi,-n = ——v-n {n^O), (56) 
1 — p n 

«i = 77o, (57) 

Wifi is free. (58) 

Nextly, by setting k = 2 (1551) gives us rjw^ = u>2 — u?2 + /«iu>i + Therefore we have 

1 / v l <r l 



W 2 ,-n = 1 _ ( ^ -j—plV-n+lV-l + V-nWl,0 ~ «1 ^ _ ^-n ) 7^ 0), (59) 



^— ' 1 — 7/ 



7^2,0 is free, (61) 
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and so on. In this way, we obtain a series of constraints (!52|) - (l54|) etc, besides the relations 
giving the nonzero Fourier modes of w^s in terms of the ?7fc's. 

Conjecture 3.12 For k = 1, 2, ■ • • and nonzero integer n, Wk- n is expressed in terms of 
the r/k 's and u^o, • ■ ■ , Wk-ifl- 

Now we present a conjecture that the constraints fl52|) - fl5^|) etc can be written as sums 
of multiple integrals whose kernel is simply given by products of theta functions with the 
period 5: 

= -i p V* e ™J2(-l) n p nin+1)/2 e 27rinx (62) 

neZ 

= ip 1/8 e- mx JJ(1 -p n - l e 2 ™ x ){l -p n e~ 2mx ){l -p n ). 

n>\ 

Note that we have ■&[ = i?i(0) = n„>i( 1 -p n ) 3 and -&i(-x) = -$(x). 

Definition 3.13 We define the quantities I n (n = 1,2, • • ■ ) by the multiple integral 



, , dxi ■ ■ ■ dx n x (63) 

x n T7 a X ~f- ?-vM---vM, 

i<r<7<n ~ x i + - Xj - V 

where we have used the theta function with the period 5 and denoted $i(x) = $i(x; 5) for 
short. 

Remark 3.14 An explanation is in order here. Similar multiple integrals as in (E3j) were 
found in J7|/ when one of the authors studied the family of Macdonald difference operators 
in terms of the free field construction. In JE/, it is found that the algebra found by Feigin 
and Odesskii |17j/ is the underlying algebraic structure for such integrals, including some 
elliptic extension of the Macdonald theory. From this point of view, one can regard the 
integrals If63\) as a certain classical limit of this. Some detail will be given in Section^ 

Conjecture 3.15 Let the degree of be k. The I n 's are homogeneous polynomials of 
degree n in the K k 7 s. Namely the quantities I n ? s are conserved under the time evolution 
described by the integro- differential equation (Tjp. 

Examples: 

h = Ki, (64) 

h = K 2 + -n « ls (65) 

Vn^O V / 

h = «3 + X _ p n K ^ + ( 66 ) 



1 y Q -P q q -PQ _ 2 y q (g ~P q ) 
6\ l-p m l-p n ^ (l-p n ) 2 ' lJ 
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and so on. 



3.4 Conserved densities 



In this subsection, we continue our study on structure of the conserved quantities from 
the point of view of conserved densitiess. 

First let us study the operator WS k W~ 1 in some detail. Set 



k-l 



(67) 



Lemma 3.16 Write u = l,Uko = 1 for simplicity. Then we have 



1=0 



UklW, 



(k-l) 
p-l 



w r , 



Hence Uk \ can be expressed as the determinant 



Ubl 



(fc) 
w[ ' 




1 



( fc ) „„(fc-i) 



Wo W 



(k) (k-l) 
W l W l-1 



1 

Wi 

Wl-i 



w 



(1) 



((it 



(69) 



'l wi 

As usual, we write res A = ao for any A = Y^=k a j^^ ■ 
Proposition 3.17 For any positive integer k, Uk t k = res (WS h W~ 1 ) is a total q- difference. 

Proof. For j — 1, 2, 3 - • • , set 



w 



(j) -, 
l L 

U) ..,(3-1) 











Wo w 



U) (i-i) 



w 



(i) 



and set £q = 1. With this notation we can write the determinant ( 1691) written for I = k as 



Uk,k 



(-i) fe E(- 1 )'(^-- w N i) 



3=1 



□ 
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Here are some examples of u^f. 



«1,1 


= res^WSW' 1 ) = (1 - 5)^, 


(70) 


Mo i 


= (1 - S 2 W, 


(71) 


«2,2 


= resiWS^- 1 ) = (I- S)(w 2 + w 2 -w 1 w 1 ), 


(72) 


«3,1 


= (1 - S>i, 


(73) 


«3,2 


= W 2 - W 2 + WiWi - WiWx, 


(74) 


«3,3 


= TesiWS 3 ]^- 1 ) = 


(75) 




= (1 - S)(w 3 + W 3 + W 3 - WiW 2 - WiW 2 - WiW 2 + WiWiWt), 





etc. 

Now we can state the structure of the residue of Bk = {W(f)(S) k W~~ 1 ) + . 

Proposition 3.18 For any positive integer n, there exists a difference polynomial of 
the Wk 's such that resB n = const. + (1 — S)E n . 

Examples: 

res Si = 1 + (1-S)w 1 , (76) 
res5 2 = (0 2 + 20 2 ) + 20 1 (l-S)w 1 + (l-S)(w 2 + W 2 -u> 1 wi), (77) 
resS 3 = {4>\ + 60i0 2 + 30 3 ) + 3(0? + 2 )(1 - S)w 1 + 30i(l - S)(w 2 + w 2 - w x w x ) + 
+ (1 - S)(w 3 + w 3 + w 3 - w 1 w 2 - WiW 2 - wi%v 2 + W1W1W1), (78) 

and so on. 

Now we assume the reduction condition (13 



Lemma 3.19 From the evolution equation d tm L = [B m ,L], we have 

d tm r ] = r ] -(l-TS- 1 )resB m . (79) 

Hence from Proposition \3.18[ we have 

d tmV = r ] -(l-TS- l )(l-S)E m . (80) 

Suppose that the £ m 's are holomorphic on D and periodic E m (x + 1) = E m . Then 
from Corollary 12.21 and (18"U|) we have some quantity H m which satisfies d tm r] = n ■ T(H m ). 
We have a conjecture which explicitly gives H m . 

Definition 3.20 We define the densities H n (x) (n = 1, 2, ■ ■ • ) by the n — 1-fold multiple 
integral 

Hn{xi) = (_i)n-i_!_ ... [ 1/2 dx2 ... d x n x (81) 

x J] T7 ~t~^\ q X ~r ;-v(xi)---vM, 

i<r<7<n ^ \ Xi ~ x i + vM x i - x i - 7) 

namely we have nl n = J^y 2 dxH n (x). 
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Conjecture 3.21 By suitably choosing the constants 4>k's and c, we have the equality 
H m = E m — E m + c, namely we have 

d tmV = r ] -T(H m ). (82) 

We show some examples. Firstly, we have E\ — Ui, and H\ = wi — wi + c. By choosing 
the constant as c = K\ we have Hi = rj. Nextly, we have E2 = W2 + w 2 — uiiWi + 2<f) 1 wi. 
Thus choosing the constants as 

c = 2k 2 — k>i(ki — 20i), (83) 

we have 

E 2 - E 2 + c = — (wi - wi + Ki)(«Ji - % + Ki — 20i) = H 2 . (85) 



4 Poisson structure 

In this section, we study a Poisson structure which is obtained as a certain deformation 
of the algebra associated with the family of Macdonald difference operators D r n (q,t) (see 
(3.4) r in VI. 3 of [6]), and study the relations with the integro-differential equation ([1]) and 
the Lax formulation which we have developed in the previous section. 

One of the authors studied [7] the family of Macdonald operators in terms of a Heisen- 
berg algebra and its Fock representation (namely in the case of infinitely many variables 
x\, x 2 , ■ ■ ■)■ In [8] is given a systematic description of the relation between this and the al- 
gebra obtained by Feigin and Odesskii [T7] (whose trigonometric limit to be very precise). 
In other words, we analyzed the commuting family of operators introduced by Macdonald 
from the point of view of a pairing between the two kinds of quantum i.e. noncommutative 
algebras. 

Once we note that the algebra of Feigin and Odesskii is originally constructed over 
an elliptic curve (which contains three parameters in this most general setting) one may 
easily find the corresponding elliptic deformation of the commutative family which are 
acting on the Fock space containing three parameters, say q, t and p. 

It is interesting to note that these underlying algebras, the Heisenberg algebra and 
the Feigin- Odesskii algebra, become commutative in the limit t — > 1 with q and p fixed. 
Hence by setting t = e h and considering the limit H — > 0, one may naturally define Poisson 
algebras on the corresponding commutative algebras. 

Because we lack the space, we skip the derivation and just give the resulting Poisson 
algebra, then compare the relations with those coming from the Lax formalizm. We again 
use the notations q = e 2mi ,p = e 2nlS for simplicity of display. The Poisson algebra we 
study is generated by {\ n \n £ Z \ {0}} with the Poisson bracket 

(1 - q n )(l -p n q~ n ) 

{\) ^m} = : ~ fin+m,0, (86) 
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where 5 m n denotes the Kronecker delta. Let e be a constant and set 



n{x) = J2v n e~ 2mnx = eexp ( ^ \ n e~ 2mnx j . (87) 

neZ Vn^O / 

Proposition 4.1 We have 

Let us define I n and H n (x) by the same equation (1631) in Definition 13.131 and ( 1HTT) in 
Definition 13.201 respectively from the quantity defined by flHTj) . Here is a crucial remark: 
we are not working with the dependent variables described by the Lax formalism, but 
we are starting from the Poisson algebra and trying to reconstruct the same hierarchy 
described by the Lax operator L together with the reduction condition (EMI) . First, one 
can prove the following. 

Proposition 4.2 For any positive integers k and I, we have 

{h,li} = 0. (89) 

Let I n be our n-th Hamiltonian, and set dn/dt n = {I n ,n(x)}. 
Proposition 4.3 We have dn/dt n = {I n ,n(x)} = n{x)T(H n ). 

For example, we have 

—n{x) = { V o, n(x)} = r,(x) £ * — q -f- / q V ne 2 — ■ (90) 

The RHS is nothing but r/(T^), hence we recover the integro-differential equation ([I]). 
Note that the time evolution in general takes the same form as we conjecture for the Lax 
formulation (see Conjecture I3.2ip . 

We note that the r also can be presented as a kind of vertex operator. Set 



r{x) = exp 



S (1 _ q n)(\ - rq - n) X n e 2nm j ■ ( 91 ) 



Then from (18711 we have 



We have 



T}(x)=e=^. (92) 

T T 



W4 r(y)} = -J2 Y^ e ~ 27riniX ~ V) I&WV)- ( 93 ) 
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From this we have 

{r]{x), r{y + S)}/r{y + 5)- { V (x), T{y)}/r{y) = (5(y - x) - 1) r]{x), (94) 
where we have used the notation 5(x) = J2 n e2mnx - This gives us 

D\T ■ t = err — r] rr, (95) 
which is exactly the bilinear equation (jEJ). 

A Weierstrass ( function 

The Weierstrass £ function ((z) = ((z; 2u>i, 2u>2) is defined by 

1 1 z 



to n 



(96) 



z — 2mwi — 2nuj 2 2muj\ + 2nu2 {2muj\ + 2nui2) 2 
where the symbol Y2 m n means the summation over (m,n) G Z 2 \ {(0,0)}. We have 

a-z) = -c(z), 

({z + 2lj 2 )=({z) + 2((u J2 ), 

and the Fourier expansion 

C(«) = M + 7^— cot 7^— + — >.i ^ sm > 97 

n I 

where p = e 27riW2 / W1 . 
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